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ABSTRACT 

In this study we investigate the relation between stellar mass, dust extinction and star formation 
rate (SFR) using ~ 150,000 star-forming galaxies from the SDSS DR7. We show that the relation 
between dust extinction and SFR changes with stellar mass. For galaxies at the same stellar mass dust 
extinction is anii-correlated with the SFR at stellar masses < 1O 1O M0. There is a sharp transition 
in the relation at a stellar mass of lO lo Af . At larger stellar masses dust extinction is positively 
correlated with the SFR for galaxies at the same stellar mass. The observed relation between stellar 
mass, dust extinction and SFR presented in this study helps to confirm similar trends observed in 
the relation between stellar mass, metallicity and SFR. The relation reported in this study provides 
important new constraints on the physical processes governing the chemical evolution of galaxies. The 
correlation between SFR and dust extinction for galaxies with stellar masses > 10 10 Mq is shown to 
extend to the population of quiescent galaxies suggesting that the physical processes responsible for 
the observed relation between stellar mass, dust extinction and SFR may be related to the processes 
leading to the shut down of star formation in galaxies. 

Subject headings: ISM: abundances — ISM: dust, extinction — galaxies: abundances — galaxies: ISM 
— galaxies: star formation 



1. INTRODUCTION 

Dust is a fundamental constituent of galaxies. Dust is 
formed from processed stellar material returned back to 
the interstellar medium (ISM) either through supernovae 
or stellar winds. Massive stars (> 8M Q ) which end their 
lives as Type II supernovae (SNe) and the AGB phase 
of intermediate mass stars (1 < Mq < 8) are consid- 
ered to dominate stellar dust production in star-forming 
galaxies while dust in the ISM may also be formed in 
situ from ac cretion of enriched gas processed by stars 
(Dwek 1998). Dust is formed from metals and therefore, 
not surprisingly, a strong correlation is observed between 
dust and the gas-phase oxygen abundance both in the lo- 
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cal universe dHcckma n et al.l Il998t i Boissier et al. 
Asari et al.1 l2007t i Garn fc BestJ l2010t Kiao et alJ 
Zahid et al.1 12012a ) and at high redshifts (jReddv et al.1 
20101). 



Leaueux et al.1 
tween stellar mass 



(|1979l ) first observed 
and gas-phase 



a relation be- 
oxygen abun- 



dance in star-forming g al axies. Using ~ 53,000 
galaxies, iTremonti et al.l (|2004l ) have since estab- 
lished a tight relation between stellar mass and the 
gas-phase oxygen abundance for star-forming galax- 
ies in the local universe. This so-called mass- 
metallicity (MZ) relation has b e en observed at lo w 
stellar masses (|Lee et al.1 120061: iZahid et al.l [2012 
and out 
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to high redshifts dSavaglio et al.l 

Erb et al.1 120061 : iCowie fc Bargerl 120081: iMaiolino et al 
~T2009I: 



2008; Mannucci et al.l 
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Zahid et al.l 120111: iMoustakas et al.l 120111: lYabe et ail 
2012f) . The metallicity is strongly correlated with stellar 



mass and the shape of the MZ relation is relatively con- 
stant with redshift. Over cosmic time the metallicities of 



galaxies at a fixed stellar mass evolve as galaxies become 
more enriched at late times. 

The MZ relation is shaped by several important phys- 
ical processes. Oxygen, the most abundant metal in 
the ISM, is primarily produced in massive stars which 
end their lives as Type II SNe, subsequently returning 
enriched material back to the ISM. However, the ob- 
served gas-phase oxygen abundance is also subject to 
large scale gas flows. Pristine inflowing gas and enriched 
outflows can both reduce the gas-phase abundance within 
a galaxy. In pristine inflows the metal content is diluted 
whereas outflows physically remove metals from the ISM. 
If the outflowing gas is enriched to levels beyond the am- 
bient ISM, either due to the direct escape of metal-rich 
ejecta from SNe or to preferential entrainment of met- 
als in ga lactic winds, the avera ge galaxy metallicity will 
decline. ITremonti et al.l (|2004f ) argue that the enriched 
outflows which more easily escape the shallow potential 
wells of low mass galaxies are responsible for the observed 
MZ relation. Because both inflows and outflows have a 
similar observational consequence, it has proven difficult 
to disentangle the effects of gas flows from observations 
of metallicity alone (see iDalcantonl 12007). In this con- 
text, the dust content of galaxies may provide important 
leverage in breaking the degeneracy of these two effects 
owing to the fact that the observed extinction in galaxies 
is dependent only on the amount of dust along the line-of- 
sight and cannot be diluted by inflows of pristine gas. In 
this study we examine the relation between stellar mass, 
dust extinction and SFR for star-forming galaxies in the 
local universe in order to better constrain the physical 
processes responsible for chemical evolution of galaxies. 

In addition to the MZ relation, a tight relation between 
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the stellar mass a nd SFR of galaxies i s observed to exist 
out to z ~ 2.5 jNoeske et al.|[2007t lElbaz et all 120071; 
Daddi et all 120071: iPannella et all l2009t iWhitaker et al.l 
20121 ) . The slope and scatter of the stellar mass-SFR 
(MS) relation are constant and independent of redshift 
and the overall normalization evolves such that at a fixed 
stellar mass galaxies at later times have lower SFRs. The 
fixed slope and scatter in the MS relation suggest that 
quiescent processes such as cosmological gas accretion 
are largely responsible for stellar mass growth since z ~ 
2.5. Understanding how the scatter in the MS relation 
is populated and how quiescent galaxies move off the 
MS relation will provide important constraints for galaxy 
evolution. In this study we examine the dust properties 
of galaxies along the MS relation to shed light on this 
issue. 

The MZ relation and its second parameter dependen- 
cies have been investigated by several groups. Most 
notable is th e relation be t ween stellar mass, metallic- 
ity and SFR. lEllison et all (120081) show that there exists 
a correlation between metallicity and specific star for- 
mation rate (sSFR) for galaxies at a fixed stellar mass. 
The relationship between stellar mass, metallicity and 
SFR w as subsequently investigated by iMannucci et al.l 
(2010) who show that at a fixed stellar mass the SFR 
is cmii-correlated with metallicity. They argue for a 
"fundamental metallicity relation" between stellar mass, 
metallicity and SFR. The lower metallicities observed 
in star-forming galaxies at early times are balanced by 
the higher SFRs in these galaxies such that the "fun- 
damental metallicity relation" do es not evolve out to 
z ~ 2. lLara-Lopez et all (|2010[) have independently 
found a "fundamental plane" relating the stellar mass, 
metallicity and SFR of galaxies which appears to match 
the observational data to z ~ 3.5. 

The observed relation between stellar mass, metallicity 
and SF R does depend on methodology and sample se- 
lection. lYates et al.l (|2012f) reexamine the "fundamental 
metallicity relation" and find that while at lower stellar 
masses the SFR is cmta-correlated to metallicity, the rela- 
tion reverses at higher stellar masses such that a positive 
correlation is observed. I Yates et al. (2012) argue that the 
"twist" in the relation is a result of gas-rich mergers at 
higher stellar masses which fuel a starburst leading to gas 
exhaustion and quenching of star formation. Subsequent 
gas accretion at levels too low to efficiently form large 
amounts of stars leads to metallicity dilution in these 
systems thus giving rise to the observed correlation. 

In order to shed light on the stellar mass, metallicity 
and SFR relation and to understand the physical prop- 
erties of galaxies populating the MS relation we examine 
the relation between stellar mass, dust extinction and 
SFR. In Section 2 we describe our sample and in Section 
3 we present our results. We provide a detailed discussion 
of selection and aperture effects in Section 4. In Section 5 
we provide a brief discussion and in Section 6 we summa- 
rize the main results of the paper. Throughout this work 
we adopt the standard cosmol ogy (.Hp, f E m , O a) = (70 
km s- 1 Mpc" 1 , 0.3, 0.7) and a lChabrierl (l2TjQl IMF. 

2. DATA AND METHODS 

We draw our sample from the SDSS DR7 which con- 
sists of ~ 870, 000 un ique galaxies spanning a redshift 
range of < z < 0.7 (|Abazaiian et al.| [2009). The sur- 



vey has a Petrosian limiting magnitude of rp = 17.8 cov- 
ering 8,200 deg 2 . The spectra have a nominal spectral 
range of 3900 - 9100A and a spectral resolution of R ~ 
2000. Both the stellar masse s, which are determined from 
the ugriz-band photometry (Stough ton et al.ll2002f) . and 
the emission line fluxes are measured by the MPA-JHU 
groupS. We adopt the DR7 values in this work but sub- 
tract 0.2 dex from the stellar masses for consistency with 
our previous work where stellar m asses are estimate d 
using a different set of routines (see Zahid et alJl20lH) . 
The SFRs in the DR7 ar e derived using the technique of 
iBrinchmann et al"1 (120041) w ith additional improvements 
given bv lSalim et al.l (|2007[ ). The SFRs are determined 
from fitting prominent emission lines in the spectra with 
the largest contribution coming from Ha and H/3 and are 
corrected for dust and aperture effects. 

The emission line fluxes are measu red by the 
MPA/JHU group (see lTremonti et al.ll2004T ). Balmer ab- 
sorption is prominent in the atmospheres of A stars. For 
integrated spectra of galaxies a correction for Balmer ab- 
sorption is required in order to estimate the emission 
line strength of Balmer lines. The emission lines are 
continuum subtracted and corrected for stellar absorp- 
tion by fitting a linear combination of the Chariot & 
Bruzual 2008 stellar population synthesis models (Char- 
lot & Bruzual, in prep) . We have scaled the emission line 
uncertainties of Ha and H/3 by 2.473 and 1.882, respec- 
tively, as recommended by the MPA/JHU group. 

From the parent sample, we select a pure star- 
forming sample of emission line galaxies. We first dis- 
tinguish star-forming galaxies from AGN by constrain- 
ing the ionizing radiation source using the [OIII1A5007, 
JNII1A6584, H/3 and Ha em i ssion lines ([Baldwin et all 
[19811: iKauffmann et al.l 120031: iKewlev et all 120061) . Fol- 
lowing IKewlev et all (|2006n . we remove galaxies where 

log([OIII]/H/3) > 0.61/ (log([NII]/Ha) - 0.05) + 1.3. 

(1) 

This selection yields a sample of 388, 000 galaxies. In or- 
der to obtain a robust estimate of the Balmer decrement 
we require that the signal-to- noise (S/N) of the Ha and 
H/3 line be greater than 8. These selection criteria yield 
a sample of ~ 157, 000 star-forming galaxies. We refer 
t o this selected samp le as the SN8 sample. 

I Groves et al.l (|2012l ) find that the H/3 equivalent widths 
and line fluxes may be systematically underestimated due 
to an overcorrection for H/3 absorption. They argue that 
this systematic underestimate in the H/3 line flux leads 
to a 0.1 mag overestimate of A v . In this study we are 
investigating the relation between ste llar mass, dust ex- 
tinction and SFR. iGroves et all (|2012f l conclude that the 
Ha and H7 lines do not suffer from the same system- 
atic errors in the absorption correction. In order to asses 
whether the effect of this possible systematic uncertainty 
qualitatively changes the relation between stellar mass, 
dust extinction and SFR, we also determine dust extinc- 
tion from the Ha/H7 ratio. We require a S/N> 3 in the 
H7 line when determining the dust extinction from the 
Ha/H7 ratio. Most galaxies in the SN8 sample (> 99%) 
have a S/N > 3 in H7. We determine the dust extinction 
from the Ha/H7 ratio in Section 4.1 in order to assess 
any systematic effects due to improper subtraction of H/3 

1 http:/ /www. mpa-garching.mpg.de/SDSS/DR7/ 
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Fig. 1. — The distribution of a) stellar mass, b) SFR (M yr 1 ), c) 
sample. 

Balmer absorption. 

We measure dust extinction from the Balmer decre- 
ment. For case B recombination with electron temper- 
ature T e = 10 4 K and electron density n e — 10 2 cm~ 3 , 
the intrinsic Ha/H/3 and Ha/ H7 ratio are expected t o 
be 2.86 and 6.11, respectively (jHummer fc Storevlll98l . 
We obtain the intrinsic color excess, E(B— V), and the 
co rrection for du s t atte nuation using the extinction law 
of lCardelli etHI (p89h and a corresponding R v = 3.1. 
We note that the results of this study are largely indepen- 
dent of our choice of extinction law as the relation pre- 
sented only relies on relative values of extinction. From 
the color excess we determine the visual extinction mea- 
sured in magnitudes from A v — R v E(B— V). 

In Figure [T] we show the distribution of stellar masses, 
SFRs, Balmer decrement and fiber covering fractions for 
the SN8 sample. We estimate the g-band fiber cover- 
ing fraction, / c , by comparing the photometric and fiber 
g-band magnitude. The g-band covering fraction is an es- 
timate of the fraction of the galaxy luminosity contained 
within the fiber. The median covering fraction for the 
SN8 sample is 0.24 (see Figure [T]i). We determine that 
selection and aperture effects do not significantly bias 
the observed relation between stellar mass, dust extinc- 
tion and SFR presented below. In Section 4 we discuss 
selection and aperture effects in detail. 

3. THE DUST AND METALLICITY PROPERTIES OF 
STAR-FORMING GALAXIES 

In Section 3.1 we present the relation between stel- 
lar mass, dust extinction and SFR. For comparison, we 
present the relation between stellar mass, metallicity and 
SFR in Section 3.2. 



3.1. The Stellar Mass, Dust Extinction and SFR 
Relation 

In Figure [5] we show the relation between stellar mass, 
dust extinction and SFR. Hereafter, we refer to the stel- 
lar mass, dust extinction and SFR relation as the MDSR. 
The data are first sorted into 16 equally populated bins 
of stellar mass and then each mass bin is sorted into 11 
equally populated bins of SFR. Each bin contains ~ 890 
galaxies. In Figure [2^i the different color curves corre- 



Balmer decrement and d) g-band fiber covering fraction for the SN8 

spond to undecilet0 of the SFR as a function of stellar 
mass. In Figure [2b we show the median dust extinc- 
tion sorted into bins of stellar mass and SFR. Again the 
curves are color coded to match the undeciles of SFR 
shown in Figure [2^ (e.g. red curves correspond to the 
median dust extinction in the highest SFR bin and the 
black curve the median dust extinction in the lowest SFR 
bin in each stellar mass bin, respectively). The median 
la scatter of the Balmer decrement within each bin is 
0.42 with 0.32 attributable to observational uncertainty. 
Given the large number of data points within each bin, 
the median standard error of the Balmer decrement for 
each bin is 0.01. 

There are several notable features present in Figure 
[21 There is a general trend for the extinction to in - 
crease with stellar mass (e.g. iBrinchmann et aLl [2004) . 
Older, higher stellar mass galaxies typically have greater 
extinction which is most likely due to the greater num- 
ber of stars in these galaxies evolving through the AGB 
phase or ending their lives as supernovae (see IDwekl 
1998). Perhaps more interesting is the relation between 
extinction and SFR at a fixed stellar mass. At stellar 
masses < 10 10 M Q the dust content of galaxies is anti- 
correlated with the SFR such that galaxies with high 
SFRs tend to have less dust extinction. At a stellar mass 
of ~ 1O 1O M0 (or A v ~ 1.2) there is a sharp transition and 
at larger stellar masses the extinction is positively cor- 
related with the SFR. Figure shows two projections of 
the 3-dimensional MDSR. The reversal of the trend be- 
tween SFR and dust extinction at a fixed stellar mass can 
be thought of as a twist in the two dimensional surface 
defining the relation between stellar mass, dust extinc- 
tion and SFR. 

The twist is also present in the data without binning 
in SFR. In Figure [3] we plot the correlation coefficient 
between SFR and dust extinction in bins of stellar mass. 
The data are binned into 16 equally populated bins of 
stellar mass, same as in Figure [2] In the stellar mass 
range of 8.5 > log(M*/M Q ) > 10 the data show a neg- 
ative correlation. At 10 10 M Q there is a transition to 
a positive correlation. We demonstrate this visually for 

2 Each of eleven equal groups into which a population can be 
divided according to the distribution of values of a particular vari- 
able. 
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Fig. 2. — The observed relation between stellar mass, dust extinction and SFR (M@ yr -1 ). a) Undcciles of the SFR as a function of 
stellar mass, b) The median Balmer decrement and visual extinction (in magnitudes, see text for details) sorted into bins of stellar mass 
and SFR. The colors correspond to undeciles of the SFR shown in a). The black error bars show the median lcr dispersion of the data in 
each bin and the red error bars show the observational uncertainty. 
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Fig. 3. — a) The Spearman rank correlation coefficient between SFR (Mq yr _1 ) and dust extinction in 16 equally populated bins of 
stellar mass. The dotted line marks the zero point, b-d) Dust extinction plotted as a function of SFR for galaxies in three of the stellar 
mass bins shown in a) . The range of stellar masses are shown in the text of each panel. The red curves are the median dust extinction in 
15 bins of SFR and the dashed curves are the 68% contours of the data. 



the unbinned data in Figure [3)>d by showing the relation 
between dust extinction and SFR in three of the stellar 
mass bins. At the lowest stellar masses there is evidence 
for another transition to a positive correlation but the 
sparsity of data at low stellar masses does not allow us 
to dra w any strong c o nclus ions. 

Both lGarn fc Best] (|2010D and lXiao et"all (|2012f) have 
studied the dust properties of star-forming galaxies in 
the SDSS. Neither of these studies report the observed 
twist in the relation between stellar mass, dust extinction 
and SFR. However, these studies do not examine the re- 
lation between dust extinction and SFR at a fixed stellar 
mass. Both studies employ a principal component analy- 
sis (PCA) on the full sample. PCA is a useful statistical 
technique for assessing the relative contribution of var- 
ious correlated parameters. PCA finds the orthogonal 
linear combinations of the variables that maximize the 
variance in the data. However, traditional PCA tech- 



niques assume a linear relation between the variables. 
The twist is not obvious using PCA because the relation 
of dust extinction and SFR is opposite at lower stellar 
masses as compared to higher stellar masses, thus "can- 
celing out" in PCA. PCA performed on a sub-sample of 
the data in restricted range of stellar mass does reveal 
the twist between dust extinction and SFR (see Figure 

El. 

3.2. The Stellar Mass, Metallicity and SFR Relation 

In FigureH^, we show the relation between stellar mass, 
metallicity and SFR. Hereafter we refer to this relation as 
the MZSR. The MZS R is determined u sing the sample 
and methodology of lYates et al.l (|2012f ). The data are 
binned into mass bins of width 0.15 dex and SFR bins of 
width 0.3 dex. The m ean metallicity determ ined from the 
Bayesian method of iTremonti et all (|2004f ) are plotted 
for each bin and the different color curves correspond to 
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Fig. 4. — The observed relation between a) stellar mass, metal- 
licity and SFR (Mq yr" 1 , c.f. Figure 1 of lYates et"aT]|20TD and 
b) st ellar mass, Balmer decrement and SFR using the T2 sample 
from lYates et al.l 1(2013 1. The data are the a) mean metallicitics 
and b) mean Balmer decrements in constant width bins of stellar 
mass and SFR. The curves are color-coded corresponding to the 
different SFR bins shown in the legend (the value given for the 
SFR is the bin center). 

the various SFR bins. The center of each SFR bin is 
given in the legen d of the figure. We refer the reader to 
lYates et al.l (p012T ) for more details on methodology and 
sample selection. 

We present Figure |4] to draw attention to the qualita- 
tive similarities in the observed MZSR as compared to 
the observed MDSR (Figure [2]). At a fixed stellar mass 
there exists an emii-correlation between the metallicity 
and SFR for galaxies with stellar masses < 10 10 M©. At 
a stellar masses > 10 10 M Q the trend reverses and a posi- 
tive correlation is observed between metallicity and SFR 
at a fixed stellar mass. As can be seen in Figur e HJd, 
the tw ist in the MDSR is also present in the lYates et al.l 
(|2012l) data. 

The anti-correlation between metallicity and SFR at 
lower stellar masses is significantly stronger than the 
anti-correlation between dust extinction and SFR. We 
determine the correlation coefficient between SFR and 
metallicity for ~ 7400 galaxies in the stellar mass range 
of 9.4 < log (M*/M Q ) < 9.5. The metallicities are taken 
from the D R7 and are determine d using the Bayesian 
technique of lTremonti et all (|2004f ) . The sample correla- 
tion coefficient between the aperture corrected SFR and 
metallicity is r — —0.41. Using the SFR determined from 
the observed Ha luminosity in the fiber the Spearman 
rank correlation coefficient between SFR and metallicity 
is r = — 0.15. The anti-correlation between stellar mass 
and metallicity is still pres ent when using Ha fiber SFRs 
(e.g. iMannucci et al.l[2010( ). however the strength of the 
correlation is diminished. 

Metallicity is a measure of oxygen relative to hydro- 
gen whereas dust extinction is dependent on the abso- 
lute number of absorbers within the line of sight. To 
first order, the observed dust extinction, unlike metal- 
licity, is independent of the gas fraction. The stronger 
correlation between SFR and metallicity as compared to 
SFR and dust extinction and much of the difference in 
the MZSR as compared to the MDSR seen in Figure [4] 
is likely due to a correlation between the gas fraction 
and SFR. Higher gas fractions may sustain higher SFRs 
while also diluting the metallicity, thus strengthening the 
anti-correlation between metallicity and SFR observed at 
stellar masses < 10 10 A/ Q . Measurements of gas masses 
in a large sample of star-forming galaxies should provide 
important insight into the relationship between metallic- 
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Fig. 5. — The observed relation between stellar mass, dust ex- 
tinction and SFR (Mq yr -1 ). Similar to Figure[2p but with dust 
extinction determined from the Hct/IPy ratio. 

ity and dust. 

4. SYSTEMATIC, SELECTION AND APERTURE EFFECTS 

In this section we investigate possible systematic is- 
sues with improper subtraction of H/3 absorption (Sec- 
tion 4.1), biases in the observed MDSR associated with 
our method of sample selection (Section 4.2) and sys- 
tematic effects of measuring global physical properties 
of galaxies from emission lines observed within a limited 
aperture (Section 4.3). We conclude that selection and 
aperture effects are not significant in our determination 
of the MDSR. 

4.1. Systematic Effects in the Hj3 Absorption Correction 

In Figure[5]we plot the MDSR with dust extinction de- 
termined from the Ha/H7 ratio. The relation presented 
in Figure [5] displays the same characteristics as Figure^). 
An anti-correlation between dust extinction and SFR at 
stellar masses < 1O 1O M0 and a positive correlation at 
higher stellar masses. We observe a ~ 0.1 magnitude 
greater extinction whe n determining A v fro m Ha/H/3 as 
compared to Ha/H7. I Groves ~etaU pOll find a simi- 
lar offset and by comparing SDSS DR7 data with DR4 
data. They attribute the difference in extinction deter- 
mined from Ha/H/3 and Ha/Hy to a systematic error in 
the subtraction of the underlying H/3 Balmer absorption. 
While systematic effects in subtracting the underlying 
H/? absorption may affect the absolute measurement of 
A v , we conclude that the observed twist in the MDSR 
is not affected. A comparison of Figure [5] with Figure 
[2Jd shows that a greater difference in A v is observed at 
higher stellar masses and SFRs. The overestimation of 
A v appears to be correlated with the stellar mass and 
SFR. 



4.2. Selection Effects 

We select star-forming galaxies from the parent sample 
using the BPT method which allows us to identify and 
remove galaxy spectra dominated by AGN emission (see 
Section 2). We require a S/N > 8 in the Ha and H/3 
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Fig. 6.— The SFR (M Q yr _1 ) distribution for the SN8 (black 
curves) and SN3 (red curves) samples. The median SFR in 30 
equally populated bins of stellar mass are shown by the filled circles 
and solid curves. The 68% and 95% contours of the distribution 
are shown by the dashed and dotted curves, respectively. 

emission line fluxes in order to obtain a robust estimate 
of the SFR and Balmer decrement from which we mea- 
sure the extinction. The strength of the Balmer lines 
scales with number of UV ionizing photons originating 
from mas sive stars and th erefore is a good indicator of 
the SFR ( Ke nnicuttl 1 19981) . We may bias our sample by 
selecting galaxies above a fixed S/N threshold because 
galaxies with low levels of star formation will have weak 
Balmer lines that may not meet our S/N requirement. 
This selection criteria could lead to a spurious MDSR if 
the bias is mass dependent. 

In order to investigate the bias introduced by our S /N 
requirement we compare the distribution of SFR as a 
function of stellar mass for the SN8 sample with a sample 
selected requiring a S /N > 3 in the Ha and H/3 emission 
lines. We refer to this sample as the SN3 samples. The 
SN3 sample consists of ~ 259, 000 galaxies and contains 
a factor of ~ 1.6 more galaxies than the SN8 sample. 
In Figure H] we plot the distribution of SFRs in 30 bins 
of stellar mass for the SN8 (black curves) and SN3 (red 
curves) samples. The median of the SFR distribution 
(solid curves) of the SN3 sample is typically ~0.1 dex 
lower than the SN8 sample except at the highest stellar 
masses where the difference is larger. 

We are interested in determining the MDSR for 
the star-form i ng se quence of galaxies as identified by 
iNoeske et all (|2007l and many others). The population 
of star-forming galaxies out to z ~ 2 is characterized 
by a near unity slope and constant scatter in the relation 
betwe en stellar mass and SFR that is independent of red- 
shift (INoeske et all 120071: lElbaz et al.l 12007) [Daddi et al.l 
120071: (Pannella et al.l l2009: Whi taker et alTl2012t) . In the 
SN8 and SN3 sample, the scatter in the SFR distribu- 
tion increases at higher stellar masses and a population 
of massive, low SFR galaxies is present in the distribu- 
tion shown in Figure [51 By examining the s ersic index of 
galaxi es on the stellar mass-SFR diagram, iWuvts et al.l 
(2011) show that this region of the diagram is populated 
by quiescent galaxies best described by de Vaucouleurs 
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Fig. 7.— The MDSR for the SN3 sample, 
profiles. Decreasing our S/N threshold slightly broadens 
and shifts the distribution of SFRs at all stellar masses 
and selects a greater number of massive, quiescent galax- 
ies. 

In Figure [7] we reexamine the stellar mass, dust ex- 
tinction and SFR relation using our SN3 sample. The 
greatest difference in the MDSR determined from the 
SN3 sample as compared to the SN8 sample is for galax- 
ies with stellar masses > 10 10 Mq and is attributable 
to the greater number of quiescent galaxies in the sam- 
ple. We conclude that our S/N selection criterion does 
not strongly select against the star-forming sequence of 
galaxies. The S/N criterion of the SN8 sample does not 
completely remove quiescent galaxies from the sample, 
though they constitute a small fraction of the sample. 
Thus, selection bias is not a significant. Comparison with 
the SN3 sample suggests that the downturn in the rela- 
tion between stellar mass and dust extinction observed 
in the highest mass, lowest SFR bin of the SN8 sample 
is a consequence of the presence of quiescent galaxies in 
the sample (see Figure [2^i) . This suggests that twist ob- 
served in the MDSR and the ani«-correlation between 
SFR and dust extinction is related to the shutting down 
of star formation in galaxies. 

We apply a higher S/N threshold in order to emulate 
selection effects that may be present in high redshift 
data. The "twist" in the MDSR is clearly observed in 
the MDSR when the S/N threshold for Ha and H/3 emis- 
sion is < 20. At a S/N> 25, there is no longer a "twist" 
in the MDSR. In the SDSS data 74% of the galaxies are 
removed when a S/N> 25 is required. At a fixed stellar 
mass this selection criterion preferentially removes low 
SFR galaxies. This demonstrates that incompleteness 
in SFR may bias the observed MDSR and care must be 
taken when investigating the MDSR in higher redshift 
samples. 

4.3. Aperture Effects 

We investigate aperture effects associated with mea- 
surements of the SFR (Section 4.2.1) and the Balmer 
decrement (Section 4.2.2). 

4.3.1. Star Formation Rate 

The SFRs used in this study are derived by 
the MPA/JHU group using the technique developed 
by iBrinchmann et all (|2004 hereafter B04). B04 
mo del the stellar conti n uum and absorption lines us- 
ing iBruzual fc Charlotl (|2003[ ) stellar population syn- 
thesis models. They model the e mission lines using 
the CLOUDY photo io nizati on code (|Ferlandl [l996h and 
iCharlot fc Longhettil (|200lD nebular emission models. 
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Fig. 8. — Relation between covering fraction, f c , and a) Ha 
SFRs and b) aperture corrected SFRs (M© yr — 1 ). The relation 
between Balmer decrement and c) Ho SFR and d) aperture cor- 
rected SFRs. The data are the SN8 sample in a limited mass range 
(9.4 < log(M„/MQ) < 9.5). The median covering fraction plotted 
by the solid red curve in 15 bins of a) Ha SFRs and b) aperture 
corrected SFRs. The median Balmer decrement plotted by the 
solid red curve in 15 bins of c) Ha SFRs and d) aperture corrected 
SFRs. The 68% contours of the distributions are shown by the 
dashed red curves. 

Dust attenuation is largely constrained using the Balmer 
decrement with small contributions from other emission 
lines. This procedure gives the SFR within the 3" fiber 
aperture. B04 apply a fiber aperture correction in order 
to obtain the total SFR. The correction is derived by cal- 
culating the dependency of SFR on color within the fiber. 
The SFR outside the fiber is accounted for by assuming 
that the color dependency of the SFR is the same inside 
a nd outside the fiber (see B04 for more details). 

ISalim et al.1 (|2007| ) der ive dust corrected SFRs for ~ 
50, 000 galaxies in the local universe by fitting the UV 
and optical spectral energy distribution (SED) with a li- 
brary of stellar population synthesis models. The SFRs 
determi ned from th e SED are not subject to aperture 
effects. ISalim et"ail (|2007D show that the aperture and 
dust corrected SFRs derived by B04 for the sample of 
star-forming galaxies (i.e. those with Ha detected with 
S/N > 3) agree with those determined from the UV 
and optical SED. There is a systematic difference of 0.02 
dex (when 3c outliers are excluded) and the scatter in 
the two meth ods is accou n ted fo r by the uncertainty of 
each method. ISalim et al.1 (|2007l ) conclude that for star- 
forming galaxies the UV-based SFRs agree remarkably 
well with the aperture corrected SFRs derived by B04. 
For the star-forming galaxies in the SN8 sample, the 
SFRs made available in the DR7 and used in this study 
are robust against aperture bias. 

The MDSR derived in this study is dependent on how 
the SFR is measured. To demonstrate the need for an 
aperture correction, we also determine the SFR from the 
dust corrected Ha luminosit y observed in a 3" fiber aper- 
ture using the conversion of iKennicuttl ([1998). We refer 
to the SFR determined from the observed Ha luminos- 
ity as the Fiber SFR and explicitly refer to the aperture 



corrected SFRs as the DR7 SFR. The relation between 
stellar mass, dust extinction and Fiber SFR is qualita- 
tively different at lower stellar masses when compared to 
the same relation using DR7 SFRs. At a fixed stellar 
mass, the Fiber SFRs are not anii-correlated with dust 
extinction at stellar masses < 1O 1O M0. We attribute 
the different relation observed between dust extinction 
and SFR at lower stellar masses to aperture effects re- 
sulting from the use of Fiber SFRs rather than the DR7 
SFRs. In Figure [5K we show the distribution of aper- 
ture covering fraction for galaxies as a function of Fiber 
SFR in a narrow mass range (9.4 < log(M*/M ) < 9.5). 
The fiber covering fraction is strongly correlated with the 
Fiber SFR. In Figure [8J: we show that Fiber SFRs and 
Balmer decrements are weakly (positively) correlated at 
Fiber SFRs < but show a negative correlation at higher 
Fiber SFRs. 

In Figure [5Jd and d we show the distribution of covering 
fraction and Balmer decrement, respectively, plotted as 
a function of DR7 SFRs. The Balmer decrement is not 
strongly correlated with DR7 SFRs. This interval con- 
tains ~75% of data. At higher DR7 SFRs there is a pos- 
itive correlation between the SFR and covering fraction. 
The correlation between covering fraction and DR7 SFR 
at higher stellar masses arises because the sample is not 
volume limited and we have applied a fixed S/N thresh- 
old. We observe a similar relation between stellar mass, 
d ust extinction and S FR using the volume limited sample 
of lZahid et alj (|2011D which is comprised of data selected 
in a redshift range of 0.04 < z < 0.1. No correlation be- 
tween covering fraction and DR7 SFR is observed in the 
volume limited sample in the mass range investigated in 
Figure[8] However, restricting the redshift range removes 
a substantial number of low and high mass galaxies. Be- 
cause the effect of not selecting a volume limited sample 
is not significant and does not change our conclusions, we 
do not apply this additional criterion in selecting data. 

By comparing Figure[5J: and d we show that not includ- 
ing an aperture correction for SFR results in a spurious 
correlation between Balmer decrement and SFR. Figure 
[H] highlights the importance of applying an aperture cor- 
rection to SDSS SFRs when examining global trends in 
order to avoid systematic bias in relations between SFR 
and other physical properties. 

4.3.2. Balmer Decrement 

Several studies have reported variati ons in dust ex- 
tinct ion with galactoc en tric ra di us (e.g. iHolwerda et all 



2005; Bois sier et al 



[20071: iTamura et all T2 009; 
20091) . However. iKewlev et all 
for systematic variation 



iMufioz-Mateos et al.l 
420 051) find little evidence 

between nuclear and global extinction for covering 
fractions > 20%. Here we test possible biases that may 
result from aperture effects in determining the Balmer 
decrement. 

If strong negative dust gradients exist in galaxies then 
the measured Balmer decrement will be biased towards 
larger values of extinction in galaxies with small covering 
fractions. In this case, the measured Balmer decrement 
will not reflect the global dust properties. We first test 
for any bias by plotting the Balmer decrement as a func- 
tion of the aperture covering fraction. In Figure |H] we 
plot the Balmer decrement as a function of aperture cov- 
ering fraction for data in a narrow mass range. For data 
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Fig. 9. — The Balmer decrement plotted against the aperture 
covering fraction, / c . The black data points are taken from the 
SN8 sample in a limited mass range (9.4 < log(M*/M0) < 9.5). 
The median Balmer decrement is plotted by the solid red curve in 
15 bins of aperture covering fraction, f c . The 68% contours of the 
distributions are shown by the dashed red curves. 

in the stellar mass range of 9.4 < log(M*/M Q ) < 9.5, 
the Balmer decrement is not strongly correlated to the 
aperture covering fraction. 

We also apply a more global test for potential bias in 
the Balmer decrement due to aperture effects. We divide 
the SN8 sample into two equally populated subsamples. 
Sample SI is comprised of galaxies that have covering 
fraction less than the median covering fraction of the SN8 
sample (f c < 0.24). Sample S2 is the complementary 
sample with covering fraction greater than the median 
covering fraction. Each sample has —74,000 galaxies. 

We want to examine the bias in the Balmer decre- 
ment associated only with the covering fraction. While 
in small bins of stellar mass, the aperture corrected SFR 
is not strongly correlated to covering fraction (see Figure 
[B})), the SFRs and stellar masses across the whole sample 
are correlated to the covering fraction. The correlation 
of the SFR and stellar mass with covering fraction is a 
consequence of the fact that SDSS is a magnitude limited 
survey and we have selected our sample using a fixed S/N 
threshold. Therefore we must take care to remove second 
order correlations between Balmer decrement and cover- 
ing fraction resulting from correlations between covering 
fraction, SFR and stellar mass. We do this by randomly 
selecting subsets of the SI and S2 data that are matched 
to have an identical SFR and stellar mass distributions 
(in bins of 0.05 dex width). We refer to these as the Sis 
and S2s subsamples. The Sis and S2s subsamples are 
comprised of ~ 39,000 galaxies each. Figure [TUb and b 
show the stellar mass and SFR distribution for Sis and 
S2s in gray. The stellar mass and SFR distributions of 
the Sis and S2s are representative of the SN8 sample and 
are identical for the Sis and S2s sample by design. Figure 
[TOb and d show the distribution of the Balmer decrement 
and covering fraction for the Sis (red histogram) and S2s 
(blue histogram) samples, respectively. 

In Figure [TT] we show the distribution of the logarithm 
of the Balmer decrement for the Sis (red histogram) and 
S2s (blue histogram) samples. The data are nearly log- 
normally distributed. The median observational uncer- 
tainties in the Balmer decrement for the Sis and S2s 



samples are 0.34 and 0.28, respectively. The Sis sam- 
ple (f c < 0.24) has a slightly broader distribution as 
compared to the S2s sample. The greater width of the 
Sis distribution may be due to the larger observational 
uncertainties of the Sis sample. Despite the very differ- 
ent distribution in covering fraction, the distribution of 
Balmer decrement in the Sis and S2s samples are very 
similar. If there were a strong bias in the measured 
Balmer decrement due to aperture effects, we would ex- 
pect a relative shift in the distribution of the Balmer 
decrement for the Sis and S2s samples. We conclude 
that aperture effects do not significantly bias the mea- 
surement of the dust extinction in the SN8 sample. 

The Balmer decrement measured in the fiber is a flux- 
weighted average over many HII regions. We perform 
a simple calculation to test whether the similarity in 
the Sis and S2s Balmer decrement distribution is con- 
sistent with the extinction observed in nearby galaxies. 
iMufioz-Mateos et al.l ([2009) measure the radial attenua- 
tion profile for galaxies in the Spitzer Infrared Nearby 
Galaxies Survey sample. The typical (median) radial 
attenuation gradient measured is -0.023 mags/kpc. We 
convert this gradient measured in magnitudes of extinc- 
tion into a gradient in Balmer decrement assuming the 
(jCardelli et al.lll989l ) extinction curve. We adopt an ex- 
ponential profile for Ha surface brightness with a fiducial 
scale length of 4 kpc. Using these values we determine 
the flux- weighted Balmer decrement as a function of fiber 
covering fraction. We find that for the fiducial values 
of radial attenuation gradient and Ha scale length, the 
Blamer decrement is overestimated by ~ 5% when the 
covering fraction is 5% (the minimum for the SN8 sam- 
ple) and a 4% overestimate when the covering fraction 
is 20%. The median covering fraction of the Sis and 
S2s samples is 16% and 33%, respectively. The relative 
systematic error in the Balmer decrement between these 
two covering fraction is ~ 1%. The calculated values 
for the relative error do not depend strongly on scale 
length adopted but do increase if a steeper attenuation 
gradient is adopted (e.g. relative error of 5% for an at- 
tenuation gradient of -0.1 mags/kpc). We conclude that 
distribution of Balmer decrement shown in Figure [TT] is 
consistent with the typic al (shallow) gradients foun d in 
local spiral galaxies fe.g. iMuno z-Mateos et al. 2009). 

As we have shown in this section, the measured dust 
extinction and SFRs are not significantly biased by 
aperture effects. However, because aperture effects are 
present in the SDSS sample in a complicated way, we also 
test whether the observed MDSR is effected by aperture 
bias. In Figure [T3] we plot the MDSR by applying an 
increasing minimum aperture covering fraction require- 
ment. In Figure I12a -d we apply a minimum covering 
fraction of 0, 0.1, 0.2 and 0.3 in selecting data and deter- 
mine the MDSR using -157,000, -143,000, -98,000 and 
—52,000 galaxies, respectively. A very similar MDSR is 
observed in FigurelT2"a-d. We conclude that selection and 
aperture effects are not significant in our determination 
of the MDSR. 

5. DISCUSSION 

In this contribution we have presented the MDSR for 
local star-forming galaxies. The physics of dust produc- 
tion, destructi on and evolut ion are not well understood 
(for review see Drainc 2003) and it is beyond the scope 



9 




9 10 11 
loq(M,/M ) 



-1 
log(*) 



3 4 5 6 
Ha/H|8 



0.0 0.2 0.4 0.6 0.. 



Fig. 10. — The black histograms are for the SN8 sample and are the same as Figurefl] The stellar mass and SFR (Mq yr — x ) distribution 
are identical for the Sis and S2s sample and are shown by the gray histograms. The distribution of the c) Balmer decrement and d) covering 
for the Sis and S2s samples are shown by the red and blue histograms, respectively. 
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Fig. 11. — Distribution of the logarithm of the Balmer decrement 
for the Sis (red) and S2s (blue) samples. 

of this study to make detailed considerations of these 
processes. It is well established that dust forms from 
metals and therefore it is not surprising that a strong 
correlation betwee n dust and m e tals is observed (e.g. 
iGarn fc Bestl 120101 ). iZahid et all (j2012bl ) derive a rela- 
tion between stellar mass, metallicity and dust extinc- 
tion. They show that the dust extinction increases with 
both stellar mass and metallicity. That is, dust extinction 
increases with stellar mass and at a fixed stellar mass, the 
metallicity and dust extinction are positively correlated. 
The correlation between stellar mass, dust extinction and 
metallicity is als o evident from PC A (jGarn & Best 2010; 
IXiao et al.ll20T2l) . 

The MDSR shows similar trends as the MDZR. While 
dust extinction can be easily estimated from the Balmer 
decrement, metallicity is substantially more difficult 
physical quantity to measure. Various methods have 
been developed to determine the gas-phase oxygen abun- 
dance from the emission line properti es of star-forming 
galax ies (for review of methods see iKewlev fc E llison 
2008). The most commonly used methods rely on the- 
oretically or empirically calibrating ratios of strong col- 



lisionally excited emission lines to recombination lines. 
However, these so-called strong line methods suffer 
from various systemati c uncertainties. In particular, 
IKewlev fc Ellison! (|2008f l show that the metallicity deter- 
mined for the same sample of galaxies can differ by up to 
~0.6 dex depending on the choice of calibration. Thus 
the observed MDSR presented in this study helps to in- 
dependently establish the observed MDZR presented in 
lYates et all (|2012D . 

The observed MZ relation and its second parameter 
dependencies serve as a Rosetta stone for understand- 
ing the physical processes responsible for the chemical 
evolution of galaxies. Because dust is formed from met- 
als, we favor a common physical origin for the twist ob- 
served in both the MDSR presented in this s tudy and 
the MDZR investigated by lYates et al.l (|2012h . While 
dust may be destroyed, metals can not and therefore the 
observed MDSR can not be explained by destruction pro- 
cesses unless the similarities in the MDSR and MDZR are 
taken to be coincidental. Moreover, given that metals are 
locked up in dust, we may expect an opposite trend in the 
MDSR and MDZR if dust destruction is the responsible 
mechanism since the destruction of dust grains should 
liberate the constituent metals, thus increasing the gas- 
phase abundance. The metallicities of galaxies are set 
by a balance between star formation and gas flows. Dust 
extinction is, to first order, insensitive to the infall of 
metal-poor gas and cannot be diluted. Therefore the ob- 
served twist in the MDSR and the MDZR can not solely 
be explained as a consequence of gas dilution. We con- 
sider differential mass loss in galaxies which is related to 
the SFR a natural explanation of the observed twist in 
t he MDSR and MDZR. 

iWhitaker et al.l (|2012f ) study the relation between stel- 
lar mass and SFR out to z = 2.5. They conclude that 
the dust attenuation observed in star-forming galaxies 
increases with stellar mass. They find that for galax- 
ies with stellar masses > 10 10 A/ Q dusty, blue galaxies 
populate the upper envelope of the scatter in the MS 
relation. Conversely, red low-dust galaxies have lower 
observed SFRs and are interpreted to possibly be in the 
process of shutting down star formation. They show that 
the MS relation at 1 < z < 1.5 for red galaxies has a shal- 
lower slope than the MS relation for blue galaxies. It is 




Fig. 12. — The MDSR for the SN8 sample but with a minimum aperture covering fraction requirement of a) 0, b) 0.1, c) 0.2 and d) 0.3. 



possible that at lower stellar masses the two relations 
cross such that a similar trend of a reverse in the corre- 
lation between SFR and dust attenuation at a fixed stel- 
lar mass may be present in the high red shift population 
of sta r-forming galaxies (see Figure 4 of Iwhitakc r et al.1 
|2012[) . However, due to incompleteness no firm conclu- 
sions can be drawn regarding the correlation between 
SFR and dust attenuation at low stellar masses. 

The result s of our study are cons istent with the inter- 
pretation of iWhitaker et afl (j2012f ). In particular, de- 
creasing our S/N threshold to S/N > 3 (Section 4.2) 
results in the inclusion of a significant fraction of qui- 
es cent galaxies as ide ntified by the sersic index analysis 
of iWuvts et all (|2Qlll ) . Figure shows continuity in the 
correlation between dust extinction and SFR observed at 
stellar masses > 10 10 Mq. This suggests that the physical 
mechanism responsible for the correlation between SFR 
and dust extinction at higher stellar masses and the twist 
in the MDSR and MZSR may be related to the physical 
processes leading to the shut down of star formation and 
the migration of galaxies to the red sequence. 

6. SUMMARY 

In this study we have investigated the relation between 
stellar mass, SFR and dust extinction. We conclude that: 

• Our analysis i s con sistent with the conclusions of 
IGarn fc B est (201fJ) that the strongest correlation 
in the data is the positive correlation between stel- 
lar mass, SFR and dust extinction. 

• The relation between SFR and dust extinction at 
a fixed stellar mass is mass dependent. At a fixed 
stellar mass, an anii-correlation between the SFR 
and dust extinction is observed for galaxies with 
stellar masses < 10 10 Mq. There is a sharp transi- 
tion at a stellar mass of 10 10 M Q . In galaxies with 
larger stellar masses there is a positive correlation 
between the SFR and dust extinction at a fixed 
stellar mass. 



• The relat ion between stella r mass, metallicity and 
SFR fsee lYates et~ai]|2012j) shows the same trends 
as the relation between stellar mass, dust extinc- 
tion and SFR. Unlike metals, dust can not be di- 
luted by inflows of pristine gas. The observed stel- 
lar mass, dust extinction and SFR relation provides 
important new constraints for understanding the 
physical processes governing chemical evolution of 
galaxies. 



• Quiescent galaxies are observed to populate the 
high mass, low SFR part of the stellar mass - SFR 
diagram. When including quiescent galaxies in the 
sample we find continuity in the correlation be- 
tween dust extinction and SFR at stellar masses 
> 1O 1O M0. The physical processes responsible for 
the relation between stellar mass, dust extinction 
and SFR at the high mass end may be related to 
the physical processes leading to the shutdown of 
star formation. 

In a forthcoming paper we develop a model incorporat- 
ing momentum driven outflows as a possible explanation 
for the observed relation between stellar mass, dust ex- 
tinction and SFR presented here. 
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